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The spectrum of short-wavelength magnons in two-dimensional quantum Heisenberg antiferro- 
magnet on a square lattice is calculated in the third order in 1/S expansion. It is shown that 1/S 
series for S = 1/2 converges fast in the whole Brillouin zone except for the neighborhood of the 
point k — (7r,0), at which absolute values of the third and the second order 1/5-corrections are 
^^ ' approximately equal to each other. It is shown that the third order corrections make deeper the 

Cn ' roton-like local minimum at k = (vr, 0) improving the agreement with the recent experiments and 

^.., numerical results in the neighborhood of this point. It is suggested that 1/5* series converges slowly 

C^ ■ near k = (7r,0) also for 5* = 1 although the spectrum renormalization would be small in this case 

due to very small values of high-order 1/S corrections. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.Gb 

I. INTRODUCTION 

f-i ■ Spin-i two-dimensional (2D) Heisenberg antiferromagnet (AF) on a square lattice has been one of the most attrac- 

c/2 , tive theoretical objects in the last two decades because this model describes parent compounds of high-Tc superconduct- 
ing cupratesi. A number of theoretical approaches have been proposed to describe the spectrum of long-wavelength 
elementary excitations (magnons) in quantum square 2D AF which results agree well with each other and describe 
quantitatively existing experimental data^"-. Meantime there are some surprising recent experimental and numerical 
_JL [ findings indicating that the standard theoretical approaches do not work for short- wavelength magnons for S* ^ 1. 
(-H . Thus, a rotonlike local minimum was observed at small T in the spin-wave spectrum ek at k = (tt, 0) in a number 

Q ' of recent experiments on square spin-i 2D AFs^^— . In particular, the magnon energy at k = (ti", 0) appears to be 
O I 7(1)% smaller than that at k = (7r/2,7r/2) in Cu(DC00)2 • 4D20.^ This local minimum is purely quantum effect 
as the classical spectrum of 2D AF is flat along the magnetic Brillouin zone (BZ) boundary connecting points (tt, 0) 
,_^ , and (0, tt) (see inset in Fig. [1]). The spectrum near the point k = (7r,0) is not reproduced quantitatively within the 
^ ■ second order in 1/S expansion^i^ and phase flux RVB techniques^. In the former case the second order corrections 
^5 ■ lead to a very small difference of 1.4% between e, n^ ~ 2.35858 and e, ,„ ,„^ w 2.39199 whereas in the last case 

\^ , ■' (tt.O) (7r/2,7r/2) 

l/-\ , this difference is too large. At the same time numerical computations using series expansion around the Ising limiti^ 
and Quantum Monte-Carloii describe the roton-like minimum satisfactorily leading to values £(J^qT^ ~ 2.18(1), 

*/") ' ^rw"w2i ~ 2.385(1) and eL g) ^ 2.16, eL/2 7r/2) ~ 2.39, respectively. The origin of the local minimum has not been 
^— N ■ clarifled yet. It is considered to be a signature of the spins entanglement on neighboring sites^. 

Existence of such a strong deviation of the spectrum near k = (tt, 0) from the result obtained in the second order 
in 1/5' is quite surprising because the second order corrections are much smaller than the flrst order ones in the 
whole BZ even for S = 1/2 (see Refi^ and below) and one could expect a small contribution from high-order terms. 
Moreover, it is well known that 1/S series for staggered magnetization, transverse susceptibility, ground state energy 
$-H ' and spin-wave stiffness of 2D AF calculated up to the third order in 1/5 converge surprisingly fast even for 5 ~ 1 
^ despite the absence of a small parameter in the theory»ii^ii2rJ^ As a result the quantitative agreement is very good 
between 1/5 expansion, numerical results and experiments. It should be stressed that quantum renormalization of 
these quantities is considerable for 5^1. For instance, quantum fluctuations reduce the staggered magnetization in 
spin-i 2D AF from its bare value of 0.5 to about 0.3. Meantime this renormalization is described quantitatively by 
the flrst few terms of 1/5 series. 

We present in the present paper results of the spectrum ej^ calculation in the third order in 1/5 and demonstrate 
that 1/5 series converges very fast in the whole BZ except for the vicinity of the point k = (tt, 0) in the case of 5 '^ 1. In 
particular, we show that absolute values of the third order corrections to the spectrum at k = (tt, 0) are approximately 
equal to and only 2.5 times smaller than the second order ones for 5 — 1/2 and 5=1, respectively. Thus, our results 
demonstrate that, unlike other quantities, quantum renormalization of the spectrum near k — (tt, 0) for 5 ~ 1 is 

described by slowly converging 1/5 series. We flnd that the excitation energy in spin-i 2D AF eL p-i « 2.3241(2) is 
3.2% smaller than eL/2 n/2) ~ 2.4007(2) that improves (but still does not make perfect) the agreement with the recent 
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FIG. 1: (Color online.) Spin-wave spectrum of spin-i AF along high-symmetry paths of the Brillouin zone shown in the 



inset. Here e\^ indicate the spectrum calculated within i-th order in 1/S so that i = corresponds to the classical spectrum 
(|10|) . Results of the series expansion around the Ising limit and Quantum Monte-Carlo (QMC) computation (available only for 
k — (7r,0) and (7r/2,7r/2)) are also shown which were taken from Ref.^'^ and Ref.— , respectively. The former results describe 
quantitatively the spectrum observed experimentally^ in Cu(DC00)2 ■ 4D2O. 



experiments and numerical results (see Fig.[T]). We suggest that despite the slow convergence of 1/S series the overall 
renormalization of the spectrum for S — 1 might be small due to very small values of high-order 1/S' corrections. 

The rest of this paper is organized as follows. We present basic transformation of the Hamiltonian and describe the 
technique in Sec. HIl Spectrum renormalization is discussed in Sec. IIIII Sec. IIVI contains our conclusion. Expressions 
for self-energy parts in the third order in 1/S are presented in an appendix. 



II. BASIC TRANSFORMATIONS AND TECHNIQUE 



The Hamiltonian of the Hcisenberg AF on a square lattice with interaction between only nearest neighbor spins 
has the form 






(1) 



We put exchange constant J = 1 in all particular numerical calculations performed in the present paper. It is 
convenient to represent spins components in the local coordinate frame using Dyson-Maleev transformation in the 
following way: 



Sj = Sp + {S^y + S'z)e'''° 



knR.1 




(2) 



(3) 



Sf = S-aU 



■3"' 3^ 



where a;, ij and z are unit vectors along corresponding axes and ko — (tTjTt) is AF vector. As a result one finds that 
the Hamiltonian ([T]) acquires the form H = Eq + X]m=i 'Hm, where Eq is the classical value of the ground state energy 



3 

and Tim denote terms containing products of m operators a and a^ . Tim = for odd m and one has for even m 

(4) 



k L 

H4 = -— ^ aLi(J2+3a-2 + ^302)0304, (5) 

^6 = OC/V2 / . >^i+3+4a_ia_2a3a4a5a6, (6) 



1,2,3,4,5,6 



where Jk = 2(cosfc2: + cosfc^), _Ek = 5 Jo, i?k = •S'Jk, A^ is the number of spins in the lattice, we drop index k in 
Eqs. dS]) and ([6]) and the momentum conservation laws J2i=i ki = and J2i=i ^i = ^^^ implied in Eqs. ([5]) and ([6]), 
respectively. 

Introducing Green's functions G{k) — (ak,a][.)cj, F{k) — {au,a-]<.)ui, G[k) — (aLkia-k)w and F^{k) = {oL]^, al^) uj , 
where k = (w, k), we have two sets of Dyson equations for them one of which has the following form: 

G{k) = G(0)(fc) + G(0)(A;)S(fc)G(A:) + G(o)(A:)[Bk + Tl{k)]F^{k), 

F^k) = G^°\k)i:{k)F^k)+G^°\k)[Bu + n^k)]G{k), ^'^' 

where G^^''{k) — {uj — i?k)^^ is the bare Green's function and S(A;), S(fc), n(fc) and Il^{k) are self-energy parts. One 
obtains solving Eqs. ([7]) and similar set of equations for G(fc) and F{k) 

Gik) ^ - + ^^ + ^(^) 



F^k) = - 



V{k) _ ' 
LJ + E]^ + T.{k) 
V{k) 

+ n(fc 
?(fc) 
Sk + nt(fc) 



G(fc) = - 

B^ + n(fc) 



p(fc) 



where 



P(fc) = 0.2 _(,(«)) '_f^(fc)^ (9) 

4"' = ^^^-^^ = ^-^0^-^^, (10) 

Vt{k) = ^k(S + S) -Sk(n + nt) -a;(I]-E) -ntf + SS, (11) 

G(fc) = G{—k), S(fc) = E(— fc) and ej^ is the spin-wave spectrum in the linear spin-wave approximation (classical 
spectrum). Quantity n{k) given by Eq. (jlip describes renormalization of the spin- wave spectrum square. We find 
fl{k) within the first three orders in 1/5 in the next section calculating corresponding diagrams for self-energy parts 
shown in Fig. [2l 

It should be noted that we do not use the conventional Bogolyubov transformation in the technique described to 
diagonalize the bilinear part of the Hamiltonian ((H). As a result anomalous Green's functions F{k) and F'l'(fc) arise 
and momenta lie in the chemical BZ that is twice as large as the magnetic one. Such an approach proved to be more 
convenient as intermediate calculations turn out to be more compact while the final results are equivalent to those 
obtained using the conventional approach.—"— 

III. SPECTRUM RENORMALIZATION 

Although the spectrum renormalization within the first two orders in 1/5 is well-known, we present here the 
corresponding expressions for the sake of completeness. 



First order in 1/S 




FIG. 2: Diagrams contributing to sclf-cncrgy parts in first three orders in 1/5*. Bold lines in diagrams (b), (f) and (h) denote 
Green's functions of the first order in 1/5* (i.e., Green's functions given by Eqs. (|8]) with self-energy parts calculated in the 
first order in 1/S). The bold line in diagram (e) denotes Green's functions of the second order in 1/S. Only diagrams (d) and 
(h)-(k) lead to the spectrum dispersion along the magnetic Brillouin zone boundary. 



A. First order in 1/5 



Only one diagram of the Hartree-Fock type shown in Fig. HJa) contributes to the spectrum renormalization in the 
first order in 1/S. The result can be represented in the form 



E<-"-\k) = MA + B), 
nt('^)(fc) = MA + 2B), 



.(1) 



.(0) 



1 



2{A + B] 
2S 



.(0) 



0.158 
25 



where the following two constants are introduced: 



s4 



B = - 



N ^-^ 



SJq 



i 0.2756, 

JO) 



2e 



(0) 



-0.1966. 



(12) 
(13) 
(14) 

(15) 



(16) 



(17) 



It is seen from Eq. (|15p that renormalized spectrum remains flat on the BZ boundary in the first order in 1/S because 
Jk = for |A:^| = TT - |fc^|. We draw e^^ for S = 1/2 in Fig. [l] using Eq. ([TSll. 



B. Second order in 1/S 

Diagrams (b)-(d) shown in Fig. [2] contribute to self-energy parts in the second order in 1/5. Diagram (b) is a 
schematic representation of the correction from diagram (a) of the second order in 1/5 which arises after calculation 
of the diagram (a) with Green's functions of the first order in 1/5 (i.e., Green's functions given by Eqs. ([5]) with 
self-energy parts given by Eqs. P^ - P^ ). 



1. Diagrams (h) and (c) 

It is convenient to group expressions for diagrams of the Hartree-Fock type (b) and (c) with the result 

S(''^)(fc) = 0, (18) 

n^^^Hfc) - J^^^^^^ (19) 

ntc^)!^) . J.^-t^±^, (20) 

where A and B are given by Eqs. P^ and (IT71) . respectively. As n'''^'(fc),n''^('"^)(fc) oc Jk, these diagrams do not 
contribute to the spectrum dispersion along BZ boundary. 

2. Diagram (d) 

One obtains for corrections to self-energy parts from the diagram (d) 

^^'^W = 1^ E 7 ^ ^((^Jo(ei+62 + e3)-^ei) 

ki+k2+k3=k 4eie2e3 (^w^ _ (g^ + €2 + £3) j 

X ( —S JkJiJ2J3 + 'S' •^2^3 "25* Jo>^l-k>^2'^3 
+ 5V2 J2-kJ3^3-k + ^i-k (^'-/o - £263)) 

- (ei + £2 + £3) (5V1 J2V3 J2-k + ^JkJi Ji-k (5V2 - e2e3))) , (21) 

n(.)(fc)._L 5^ ^^ + ^^ + ^3 



4616263 ( ^2 - (61 + 62 + 63 



nt(.)(fc)^ 1 E ^^ + ^^ + ^^ 



ki+k2+k3=k 4616263 f ^2 - (61 + 62 + 63) 



ki+k2+k3=k '±eie2e3 [^ - [<^l + ^2 + t3 

X ( ^■;^'5' J1J2J3 + '2'S Jo>/2>/2-k'/3 — S' JiJ2-\^J2J3 

-5JlJ2-kJ3-k(^V2-6263)), (22) 

1 

iv2 

ki- 

X ( —■^5' J\^JlJ2J3 + 25" JQJkJlJ2-ls.J3 — S J1J2J3 
+ 2S JoJiJi-kJ2 —S Ji_]^JlJ2J3 

- {2S.hJl + SJ1J2J3 + 2SJ0J1J1-1. + S.hJ2-kJ3-k) {S^Jo - £263)) , (23) 

where we drop superscript (0) in 61 2 3 to light notations. Sums in Eqs. (|2ip -(|23 p over each momentum were calculated 
numerically by summing up L^ points in BZ with some particular values of i ranging from 20 to 200 and extrapolating 
the results to L = 00 using the formula Aoo + Ai/L + A2/LP' + . . . , as it is done in previous papers. ^'^'^^ Appropriate 
symmetry of the summands was also used. 

In accordance with previous results^s^ we obtain that this diagram leads to a very small difference of 1.4% between 
t, n^ ~ 2.35858 and t, ,„ ,„> w 2.39199. It is seen from Table IT] that the second order corrections are much smaller 

than the first order ones in the whole BZ for all S. The spectrum e\^ is presented in Fig. [T] for S — 1/2. 

C. Third order in 1/S 

One has to analyze in this order diagrams shown in Fig.[2je)-(k). Diagram (e) represents the second order correction 
from the diagram (a) which should be calculated using Eqs. P^ - (P5)) . Bold lines in diagrams (f) and (h) denote 
Green's functions of the first order in 1/5*. Expressions for self-energy parts in this order are quite complicated and 
the reader is referred to [^ for some detail of their calculation. It can be shown (see Appendix [^ that diagrams of the 



TABLE I: Expressions are presented of the spin-wave spectrum ej^ within the third order in 1/S in some representative points. 
Here ej^. is the classical spectrum given by Eq. (|10|l . The corresponding values of ej^ are also shown for S — 1/2. Notice the 
smallness of the second order l/S'-corrections as compared with the first order ones for all points and all S. In contrast, the 
absolute value of the third order correction is approximately equal to the second oder one at k = (tt, 0) for S — 1/2. 



¥■» 

(^,0) 

Stt tt 

T'4 
/TT 7r\ 

I2' 2/* 
Stt Stt 

T'T 



arbitrary 5* 



,(0) 
=k 

.(0) 
=k 

.(0) 
=k 

,(0) 
=k 

,(0) 
=k 

,(0) 
=k 

,(0) 



, 0.15795 0.02476 

1+ ^^ + 



25 



(2S)2 



, 0.15795 0.02879 

1+ ^^ + 



25 



(2S)2 



, 0.15795 0.02538 

1+ ^^ + 



25 



(2S)2 



, 0.15795 0.02134 

1+ ^^ + 



25 



(2S)2 



0.15795 0.02967 

25 (2S)2 

0.15795 0.03805 

^ 2S ^ (2S)2 

, 0.15795 0.02914 

1+ „„ + 



25* 



(25)2 



0.0033(3 ) ' 

(2S)3 
0.0042(1) ' 

(2S)3 
0.0118(1) 

(2S)3 
0.0172(1 )- 

(2S)3 
0.0065(1) - 

(2S)3 
0.0043(1) 

(2S)3 
0.0005(1 )- 

(25)3 



5=1/2 



1.2290(3) 
2.0482(2) 
2.3179(2) 
2.3241(2) 
2.3622(2) 
2.4007(2) 
1.6781(2) 



Hartrcc-Fock type presented in Fig.[2je)-(g) do not change along BZ boundary so that only diagrams (h)-(k) give rise 
to the spectrum dispersion in these directions in this order. In contrast to the second order corrections ([2T |) ~ ([23| one 
has to calculate triple sums over momenta in the third order. This procedure requires pretty much computer time. 
Then, we focus on short-wavelength magnons as their spectrum renormalization is expected to be most pronounced 
and calculate e\^ in a number of points with |k|, |k — ko| > tt/8. The results are presented in Fig. [1] (for S = 1/2) 
and in Table n 

It is seen from Table U that the third order corrections are noticeable only for S' ~ 1 and only in the vicinity of the 
point k = (tt, 0). In particular, absolute values of the third and the second order corrections in 1/S are approximately 

in spin-i 



equal to each other at k = (tt, 0) for S = 1/2. The excitation energy in spin-i 2D AF e[^L « 2.3241(2) is 3.2% 



smaller than e 



(3) 



;/ ,2 7r/2"i ~ 2.4007(2) that improves the agreement with the recent experiments and numerical results 
leaving it, however, far from being perfect. Thus, our calculations demonstarte that quantum renormalization of the 
spectrum near k = (7r,0) for S = 1/2 is described by slowly converging 1/5 series. 

It is also seen from Table |T] that at k = (tt, 0) the third order correction is only 2.5 times smaller than the second 
order one for S = 1. Thus, one can expect slow convergence of 1/S series near k — (tt, 0) also for 5 = 1. Meantime 
the overall renormalization of the spectrum would be small due to very small values of high-order 1/5 terms. 



IV. CONCLUSION 



To conclude, we calculate the spin-wave spectrum of 2D AF on a square lattice in the third order in 1/5 to examine 
the convergence of 1/5 series. Within the first two orders we recover the previous results^'* showing that the second- 
order corrections are much smaller than the first order ones in the whole BZ and for all 5 (see Table|l|. Our calculation 
of the spectrum in the next order demonstrates that the third order corrections to the spectrum are much smaller 
than the second order ones in the whole BZ except for the vicinity of the point k = (7r,0) in the case of 5 ~ 1. In 
particular, their absolute values are approximately equal at k = (7i',0) for 5 = 1/2 (see Table Hand Fig. [Ij. Thus, 
our results demonstrate that, unlike other quantities, quantum renormalization of the spectrum near k = (tt, 0) for 
5 ~ 1 is described by slowly converging 1/5 series. We find that third order corrections for the spectrum improves 
the agreement with the recent experiments and numerical results in spin-2 2D AF. We expect slow convergence of 
1/5 series near k — (tt, 0) also for 5 = 1 while the overall renormalization of the spectrum would be small in this case 
due to very small values of high-order 1/5 corrections. 
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Appendix A: Expressions for the third-order diagrams 

We present in this Appendix expressions for self-energy parts in the third order in 1/ S which originate from diagrams 
shown in Fig. [2lje)-(k). Simple codes have been written in Mathematica software to generate the majority of these 
expressions. To make it compact we present below expression for the sum of anomalous self-energy parts Il(k) +Ti'^ [k) 
rather than for n(fc) and 11^ (fc) separately because only this sum contributes to the spectrum renormalization in this 
order (see Eq. pTjl ). 

1. Diagrams (e) and (f) 

It is convenient to group contributions from diagrams of the Hartree-Fock type shown in Fig. [2][e) and (f ) . One has 
after simple calculations using expressions (jl2l) -(( 



^^^^ ^ ki+k.+k3+k4=o 32eie2e3e4 (ei + £2 + £3 + £4) 

X I — oo Ji J2^3'^4^1 "^~ '^1<^4 ( "~l0O J2 -^3^1 H~ 8*^ '^0'^2 '^1+4^4 

+85^0 Ja J4 J2+4 (ei + £4) + 2J2 J3 (3eie2e3 - S^ (SJ^ + 2J^^^) 64)) 

+ 4 {J^J^+^ei {-S^J^ + £263) - iJlJi {S'^Jlei + e^ {f JI - £163)) 

-I- J2 (45* JQJiJ2+i<^2 + 2JoJ4ei yS J3 (2J4J1+4 -I- J3J2+4) — 2J2+4e2e3) 

+ S Jq J3+4 ( J3 J2+4 (ei — £4) — JiJi+i^i) 

+J4,.h+iei (-S*^ J3 J4 J2+4 + •h+i^-i^i)))] , (Al) 



^J]^J_ y- 1 

Jo 8iV3 ^ eie2e3e4 (ei + £2 + £3 + £4) 

ki-(-k2+k3+k4=0 

X (5J3 J4 {S^JlJ2 - 252 Jo Ji J2+4 J4 + ( Ji J2 - 2 Jo J3+4 + 2J3 J4) {S\Jl - £162))) , (A2) 

where A and B are given by Eqs. ((T6)) and p7)) . respectively. Sums in Eqs. (jAll) and (jA2p were calculated as it is 
described above in Sec. HUB 21 for diagram (d) with 20 < L < 96. These sums arise after taking into account in 
Green's functions involving in the diagram (e) contributions to self-energy parts from the diagram (d). We obtain 
numerically from Eqs. dJT]) and ^^ E("-^)(fc) = ^^^|y^ and n('=/)(fc) + nt('=/)(fc) = Jk^^Jg^. 



2. Diagram (g) 

Corrections to self-energy parts from another Hartree-Fock diagram shown in Fig. [2{g) have the form 






16iV3 fr;^ , . „ eie2e3e4 (d + £2 + £3 + £4) 

ki+k2+k3+k4=0 



X (5J3 J4 (5^3 J2 - 2S^JoJiJ2+iJi + {J1J2 - 2 Jo J3+4 + 2 J3 J4) (5^2 - £162))) , (A3) 

^' ^' JoN\^^2^^^ 85£l£2£3£4(£l+£2+£3+£4) 

ki+k2+k3+k4=0 

X Ji+2( "S* (—25' Jq J1J2 + 5 Jot/1+2 — 35* JoJi J3J4 + Jo ('S' J1+3 (2J2J3 + JiJa) 

-2Ji+2£l£2) + Jl J4 (5"^ J2 Jl+2 ^3 + -/l+3£2£3)) 

+ ( Jl+2£l£2£3 + 25^2 (- J3 Jl+3£l + Jo Ji£3)) £4) ■ (A4) 



Notice the identity of sums in Eqs. (KM and (|X2]) . Sums in Eqs. (TO)) and (|X4)) were calculated with 20 < L < 96. 
We obtain numerically from Eqs. (Ml and (M) S(9)(fc) = - 0-05g2(2) ^^^ u'-a^k) + nt(s)(fc) = -A °-™j^\ 

3. Diagram (h) 

It is convenient to divide corrections from diagram shown in Fig. [21[h) into two parts: E'^''^(fc) — T,-^ '(k) + Y,^ (k), 

n('')(fc) = n^''^(A;) +n^''^(fc), and nt(^)(fc) = Ul^''\k)+Ul^'^\k), where the first terms arise after taking into account 
first order 1/S corrections to self-energy parts in numerators of Green's functions in diagram (d). As a result one has 
for them 



^^"W = ]^ 



Yj 7 ^((£i+£2+£3) 

ki+k2+k3=k 4£i£2£3 ( W^ - (£1 + £2 + £3) ) 



X Q(3A + S)5Vo Jk Ji J2 J3 - (3A + 25)^2 JoVi_k(JkJi + 2J2 J3) 

+ AJk Ji Ji_k£2£3 + {3A + B)S'^JoJiJl - (3^ + 2S)5ViVi_k J2 J3 

+3(^ + B)5Vo J2 J2-k J3 J3-k + 3(A + B)S^4Jtu -iA + B) Jo J^ke2£3) 

- w£i {2ASJpi - {2A + S)2S'Jo Ji_k J2 J3 

+ 2S{A + B)( J2 J2-k J3 J3-k + JoJLk) + A5Jk Ji J2 J3)) , (A5) 



£1 + £2 + £3 

7V2 



ki+k2+k3=k 4£l£2£3 (^^^ - (£1 + £2 + £3) j 

; (6{A + B)S^J^J2J2-k - S'^Jl ( J3 ((3A + 2B) Ji J2 + 6(A + S) Ji_k J2-k) 

-2(3A + 25) J2 Jk) + i J3 (-52 Ji J2 ((9A + 4B) J2 

|-12(A + S) J|_k + (3A + 2B) J^) + 2 (AJi J2 -H 2(A + B) Ji_k J2-k + 2AJ3 Jk) £i£2) 

- 2 Jo (2(3A + 2B)52 J2 J2 J2-k + {A + B)J^ (35^ J2 Ji_kJk - J3-k£i£2))) ■ (A6) 



Expressions for E2 (fc), IIj (/c), and 112 (^) <^an be easily obtained from Eqs. (HJ), (P^ and (P5)) taking into 

''k '^k 



account the first-order renormalization of the spectrum e:' = e:'{l + {A + B)/S) and the fact that one has to put 



£j. instead of w calculating the third-order correction to the spectrum. Sums in Eqs. (jASP and (jA6p were calculated 
with 20 < L < 200. 



4. Diagrams (i) and (j) 

Grouping expressions for diagrams shown in Fig. [2l^i) and (j) one obtains 



ki+k2+k3=k 4eie2e3 ( w^ - (ei + £2 + £3) j 
X (B5Jk Ji J2 J3 + 2BSJIjI - 2(^ + S)S'Jo Ji-k J2 Ja + 2A5J2 J2-k J3 Ja-k 
+ ;5M% ('S'^-^o - £263) j + 2AS'J2 J2-k J3^3-k (5 Jo (ei + £2 + £3) + wei) 

- (ei + £2 + £3) Q(3A + 2B)S^JiJiJ3J2-w + {A + B)J^.hJi^^ {S^J^ - £263))), (A7) 

n(-)(fc) +nt(-)(fc) = -1, Y: J1±I1±I^ -(-BS^Jl.h.hJ, 

ki+k2+k3=k 4eie2e3 (^w^ - (£1 + £2 + £3) j 

+ {3A + 2S)S'Vo JkJi J2-k J3 - SBS-Vf J2 J3 + 2{3A + 2S)S'Vo Ji Ji-k J2 

- MS-^ J^k>/i>/2 J3 - 4AS'VoVi_k J2-k J3 - (4SJk Ji^ + 2BJ1J2J3 

- 2{A + B)JoJlJl^^ + 2AJi J2_k J3-k) (5V2 - £2£3)) . (A8) 



Sums in Eqs. (|A7|) and (|A8|) were calculated with 20 < L < 200. 



5. Diagram (k) 



Expressions stemming from this diagram are very cumbersome. We present here only expression for E('^'(fc) in the 
most compact form (i.e., before integration over energies) for the particular case of the momentum k lying on BZ 
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boundary (i.e., for \kx\ — tt — \kz\ and Jk = 0) that has the form 

( — 8i^lG'5G2G'3G4 Ji_k>/2«/2-k ^ 8F4G1G3G2G5 Ji Ji_k>/4-k ~ 8i^4F5GiG3G2 Jl J4>/4-k 

— 8F2G3G5G1G4 (2Ji_i^J2 + J2-kJiJ4:-k) — 8FiF2i^5G3G4 Ji Ji_k>/2+4 

— 'S)F^G2GzGiGiJi-\^J2J2+i — 8F3G1G2G4G5J1-U.J4J2+4 — 8-F'4-F5G2G3Gi J3 J4 J2+4 

— 8F2F3G1G4G5 J3 J4J2+4 — 8F1G2G3G4G5 J3_k>/4>/2+4 — 8i^2-F4^5GiG3Ji J4_k^2+4 

— 8FiF4i^5G2G3 J2 J4-k</2+4 — 8G3G5G1G2G4 J2-k>/4-k«/2+4 

— 8G1G2G3G4G5 J2-k«/4-k^2+4 — 8F1F4F5G2G3 J2-k (Jl-kJs + J4J2+4) 

— 8i^iG3G5G2G4 (JlJi_k + •^3>^4-k>^2+4) — 4i^ii^2-F3-p4-P'5 ("^l (2>/i_k + '^^2) 
+2J4 (Ji_k J4-k + 2J3_k./2+4)) 

— AF2F3FiF5Gi (4Ji_kJ2 J3 + J2-k (Jl J3 + 2J4-kJ2+4)) 

— 8F2i^5G3G4Gi Ji_k (2 J2 J4 + ^^l-k'/3+4) 

— 4F4F5G3G1G2 (4Ji_k (Ji_k + 4) + •^2-k ( Jl J2 + 4J4_kJ3+4)) 

— 8F2F3F4F5GlJ4-k (J2-kJ2+4 + Jl J5) 

— 8F2F5G3G1G4 Jl-k (J2 J4 + ^3^5) - 16i^iF3F5G4G2 Jl_k {J2^kJ4 + J3-kJ5) 

— 8F2F5G1G3G4 (2Ji_k>/2>/4 + •^l-k'^3+4 + JlJa-kJb) 

— SF3F4F5G1G2 {2Ji_k'^2 + J2J4 + Jl-kJ4J3+4 + '/4-k {J2-\l'J4 + ^l>/3+4) + J3J4J5) 

— SF3FiiF^G2Gl (2J]^_k'^2 + '/2'/4 + «/l-k>/4-'3+4 + ('/3'^4 + ^3-k>/4-k) ^/s) 

— 8FiF3i^5G2G4 ( Ji ( J2 J4 + >/l-k«/3+4 + 2 J3J5) + 2 Ji_k (>/2-k>/4 + -/s-k^s)) 
~ 8F2G5G1G3G4 (2 J^_k"^2 + J4-\i. {JlJz+4 + •^S-k'^s)) 

— 8FiF2i^4i^5G3(2J^_l^J2_k + Jl (2Jl-k^2 + •-'^4^3+4) 
+ J4 (2J2-k^4 + J2J4-k + 2J3-k^5)) 

— AFiFiFr,G3G2 (2J2_kJ4 J2+4 + 2J4-k ( J3 J3+4 + Ji-k J5) + Ji (J2 J3 + 2J4 J5)) 

— 8F2i^5GiG3G4 ((Jl J2 + 2Ji_kJ2-k) J4-k + >/l>/3J5-k) 

— 8F1F5G3G4G2 (Jl Jl-kJ4 + J2 J3 Js-k) 

— 8F2i^3F5GiG4 (Jl Jl-kJ4-k + J2-k {'h-h-k + JsJs-k)) 

— 8F4G1G2G3G5 J2-k (J2 J4-k + J3 Js-k + Jl J2+5) 

— 8F1F5G2G3G4 ( J2-k (J2 J4 + J3 J5) + (Jl J2 + 2Ji_k J2-k) J2+5) 

— 8F1F2F3F4G5 (J2 J2-k J4 + J2 J4-k + J2-k (J3 J5 + 2Ji_kJ2+5) + Jl (Jl-k J4 + J2 J2+5)) 

— 8i^ii^3G5G2G4 (2Jf_k>/2-k + J3 J4-k J5 + Jl (2 Ji_k J2 + J5 J3+5))) , (A9) 

where the momentuin conservation laws fci + A:2 + ^3 = ^i + ^4 + ^5 = fc are imphed. The expression for J^^'''{k) for 
arbitrary momentum and those for anomafous self-energy parts are much more cumbersome than Eq. (|A9p and we 
do not present them here. Corresponding sums over momenta were calculated with 16 < L < 64. 
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